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We present a new supersymmetric approach to the Kondo lattice model in order to describe simul- 
taneously the quasiparticle excitations and the low-energy magnetic fluctuations in heavy-Fermion 
systems. This approach mixes the fermionic and the bosonic representation of the spin following 
the standard rules of superalgebra. Our results show the formation of a bosonic band within the hy- 
bridization gap reflecting the spin collective modes. The density of states at the Fermi level is strongly 
renormalized while the Fermi surface sum rule includes n c + 1 states. The dynamical susceptibility is 
made of a Fermi liquid superimposed on a localized magnetism contribution. 
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Heavy-Fermion systems are typical examples of systems with strong electron-electron correlation effects. They 
are at the origin of the large effective masses as observed by specific heat, susceptibility and de Haas van Alphen 
measurements. The formation of the heavy quasiparticles is believed to result of the Kondo effect in these 
spin-fermion systems. At the same time, a serie of magnetic properties are found superimposed on the Fermi- 
Liquid type excitations. Let us quote the evidence for two distinct contributions in the dynamic structure factor 
measured in CeCug and CeRu2Si2 by Inelastic Neutron Scattering: a q-independent quasi-elastic component 
typical of the local Kondo effect, and a strongly q-dependent inelastic peak reflecting the magnetic correlations 
due to RKKY interactions. Also the onset of a metamagnetic transition under magnetic field in CeCu6 and 
CeRu2Si2, as well as the development of small ordered moments in UPt3 are characteristic of local magnetism. 
Quite obviously, the physics of heavy- Fermions is dominated by this duality between Kondo effect and localized 
moments. 

Traditionally, the spin is described either in fermionic or bosonic representation. If the former representation 
used for instance in the l/N expansion of the Anderson or Kondo Q lattice, appears to be well adapted for the 
description of the Kondo effect, it is also clear that the bosonic representation lends itself better to the study of 
local magnetism. This constitutes our motivation to introduce a new approach to the Kondo lattice which relies 
on an original supersymmetric representation of the impurity spin 1 /2 in which the different degrees of freedom 
are represented by fermionic as well as bosonic variables. The concept of supersymmetry originally stems from 
field theory. It has been applied in condensed matter physics on two occasions: first in the problem of localization 
in disordered metals ||, and secondly in the study of the underscreened Kondo impurity ( S > 1/2) |Q. In the 
latter case, only the spin channel, which is screened by the conduction electrons (s = 1/2), is treated fermionically, 
while the residual degrees of freedom are represented bosonically. 

We propose the following supersymmetric representation of the spin S = 1/2. The two states constituting 
the basis may be written: |l/2, 1/2) = (xf\ + yb\)\0) and |l/2, -1/2) = (x/| + yb\)\0), where /t and b\ are 
respectively fermionic and bosonic creation operators, x and y are parameters controlling the weight of both 
representations and |0) represents the vacuum of particles: f a \0) — b a 0) = 0. The spin operators are then given 
by S a = ^2 aa , f^T^'fa' + ^ a T^ a ,b a > = Sj + 5", where a = (+, — , z) and r Q are Pauli matrices. S/ and St, are 
similar to Abrikosov pseudo-fermion and Schwinger boson representation of the spin respectively. 

One can easily check that the new representation satisfies the standard rules of SU(2) algebra: 15,771) are 
eigenvectors of S 2 and S z with eigenvalues S [S + 1) and m, [S + , S~] = 2S Z and [S z , S^] = ±S' ± provided that 
the following local constraint is satisfied: 

J2fiU + b% = n f +n b = l, (1) 
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which implies x 2 + y 2 = 1. In order to eliminate the unphysical states as (xfl + yb*_^j |0), we need to introduce 

a second local constraint and take: Q = J2 a a (fab-a + b_ a f a ) = 

The partition function of the Kondo lattice can be written as a path integral: 

Z = J Vc ia V fivVbivdXidKi exp - f£ dr(C + U 
+ E i A i (n /i +n 6i -l)+E i ^iQi 

with C a = £ i(7 c L d rCta + fld T f ta + bl a d T b ia (2) 
and H = J2ka £ k c l <J c kcr+J2i E o(nf t +n bt ) 
-ME 4 ( n c, +n fi +n bi ) , 

where J(> 0) is the Kondo interaction, and two time-independent Lagrange multipliers Aj and Ki are introduced 
to enforce the local constraints. 

Performing a Hubbard-Stratonovich transformation and neglecting the space and time dependence of the fields 
introduced CTj, Aj, r\i in a self-consistent saddle-point approximation, we have: 

Z = J dridrfC (a , A , rj, V*) Z{q, rf) 
Z(t],r]*) = J Vc ia Vf ia Vb ia exp[- ^ dr (jC + W) 



with C (cr , A , m*) = ex P/5 (a - ^ 



(3) 



ej — fj, cr 
H = | (To £fe - M r? 

77* e/ — /U 



£/ = £ + A • 



Note the presence of a Grassmannian coupling 77 between Cj CT and bi a , in addition to usual coupling ctq between 



c; La and /io- responsible for the Kondo effect. Ho is a supersymmetric matrix of the type ^ ^ ^ ) in wliicli 0.6 
(p,cr) are matrices consisting of commuting (anticommuting) variables. 
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FIG. 1. Sketch of energy versus wave number k for the three bands a, (3, 7 resulting of the diagonalization of supersymmetric 
Ho- Note the presence of the bosonic band 7 separating the fermionic bands a and (3. Note also the flatness of the lower band 
at k — kF- 

Next step is to diagonalize the supersymmetric Ho- The resulting spectrum of eigenenergies is schematized in 
Figure ^. The new feature brought by the supersymmetric approach is the presence of a dispersionless bosonic 
band 7 within the hybridization gap separating the two fermionic bands a and /?. 

In the scheme we propose, oo and Ao are slow variables that we determine by solving the saddle-point equations, 
while 77, 1]* are fast variables that are defined from a local and instantaneous approximation. Performing the 
functional integration of ([|) over the fermion and boson fields yields || a superdeterminant (SDet) form written 
as follows 



Z(r],r]*) = SDet(d T + H) , 

where SDet(d T + H) = H) , 

G- 1 =d r + a and D' 1 = d T + b . 

Expanding to second order in 77, 77* allows us to define the propagator G m * (k, iuj n ) associated to the Grassmann 
variable r\ and hence the closure relation for x\ = (rjrj*): 

x o = p GVOM^n) , (4) 



with Gnv'^n) = {1 _ ^^.^ 

and U° cb = — ^2 G cc (k + q, iui n + iu> v )D(q, iuj n ) 



The resolution of the saddle-point equations, keeping the number of particles conserved, leads to: 

-Vf = (£/ - m) = D exp 
2p (a$ +xl) 



2J Po 



1 = 



-VF 

(°o + 4) 

D 



(5) 
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where po = 1/2D is the bare density of states of conduction electrons. From this set of equations we find : _- , 



: £7 = 

and fj, = yp. These equations constitute an extension of the usual mean-field results of the 1/N expansion []T|,p[ to 
the supersymmetric case. A key role is played by the parameter Xq determined by the closure equation (0) that 
controls the weight of bosonic over fermionic statistics. We have studied the J-dependence of Xq and reported 
the results in Figure |^. Xq first increases when J is reduced and then decreases reaching zero at J = 0. 




FIG. 2. .//Independence of the renormalized Grassmannian coupling x\ = {rjrj*), fixing the weight of fermion and boson 
statistics. 



At zero temperature, only the lowest band a is filled. Let us then discuss the density of states p{Ep) and the 
respective weight of /, c and b in the quasiparticle a formed at the Fermi level. A naive argument might lead 
us to think that p(Ep) could be depressed when gaining new bosonic degrees of freedom. In fact nothing of the 
sort happens: we found an enhancement of the density of states unchanged from the standard 1/N expansions: 
p{E F ) 1 , {ol+xl) D 



= 1 



Po 
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•C 1. The weight of f ,c,b in the quasiparticle a at the Fermi level arc given 



by uj,,vj, and p\p\ respectively. The bosonic weight behaves as Xq as a funcion of J. Note that all of these results 
are valid as long as the bosons are not condensed. 

Largely discussed in the litterature || is the question concerning the Fermi surface sum rule: do the localized 
spins of the Kondo lattice contribute to the counting of states within the Fermi surface or do they not? Depending 
on the answer, one would expect large or small Fermi surfaces. The supersymmetric theory leads to a firm 
conclusion in favour of the former. It is easy to check that the number of states within the Fermi surface is just 
equal to n c + n\, + n/, i.e. n c + 1. The Luttinger rule includes both conduction and localized electrons, even if 
the latter described as spins in the Kondo hamiltonian are deprived of charge degrees of freedom. 

Let us next consider the dynamic spin susceptibility. It consists of two contributions, one fermionic Xff smd 
one bosonic xtb- Writing them in the cigenbasis {a., (3, 7} and neglecting the interband term in the frequency 
range of interest, one gets: %(q, w) = Xaa(q, + x 77 (q, u>), with 



a = tt/2 for Ek — k 2 /2m and v* F is the Fermi renormalized velocity. The first contribution is of Lindhard-type. 
It is characteristic of fermionic quasiparticlcs. The second contribution is bosonic. It is zero at T = since the 
7-band is empty. It increases as the temperature gets larger, with a peak of x"(^)/w at u) = 0. 

In summary, we have shown that the supersymmetry proves to be a powerful tool already at mean-field level 
to account for both the Kondo effect and the localized magnetism present in the heavy-Fermion systems. We 
show how the weight of each of the bosonic and fermionic component is controled by a Grassmannian coupling. 
Following the standard rules of superalgebra, we find an increasing weight of bosonic component when the 
Kondo coupling is reduced before decreasing and reaching zero at J = 0. The key results are the prediction of 
two contributions in the dynamical spin susceptibility, of Fermi liquid and localized magnetism-type respectively, 
the existence of an enhanced density of states at the Fermi level reflecting the heavy effective mass while the 
Fermi surface sum rule includes the total number of n c + 1 states. 
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